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I.  INTRODUCTION 


One  of  the  principle  roles  of  the  Ballistic  Research  Laboratories 
(BRL)  is  to  generate  experimental  data.  Often  projectile  velocities, 
orientations  and  behind  the  target  space-time  distributions  of  fragments 
are  required.  A common  technique  for  obtaining  such  measurements  is  the 
use  of  orthogonal  flash  x-rays  or  spark  shadowgraphs.  The  basic  concept 
is  illustrated  in  Figure  1.  Two  localized  radiation  sources,  and  S2, 

project  respective  shadows,  Fj  and  F^,  of  the  projectile  on  the  plane 

photographic  plates  oriented  at  right  angles  to  each  other.  It  is  a 


straightforward  geometry  problem  to  use  these  source  and  image  coordi- 
nates  to  obtain  the  projectile's  space  coordinates.  Let  X^,  Y^  and 

Zp  be  the  projectile  coordinates  in  a Cartesian  frame,  0,  such  as  indi- 
cated in  Figure  1.  Conceptually,  one  obtains  Y^  and  in  terms  of 

from  the  S_,  F„  coordinates,  and  similarly  one  obtains  X and  Z in 
2 2 P P 

terms  of  Y^  from  the  S^,  F1  coordinates.  Thus  one  has  two  equations  in 

two  unknowns,  X and  Y . 

P P 


Quite  often,  however,  spatial  constraints  on  an  experimental  setup 
do  not  allow  one  to  construct  orthogonal  film  planes  with  the  radiation 
sources  far  enough  away  to  sufficiently  reduce  dispersion  due  to  the 
finite  source  size.  In  such  cases  it  is  often  convenient  to  construct 
film  planes  at  other  than  right  angles.  Another  example  of  a nonrecti- 
linear  setup  can  be  found  when  the  geometry  of  the  apparatus  suggests 
the  use  of  curvilinear  coordinates.  In  many  small  arms  ranges  safety 
requirements  dictate  that  the  flight  of  the  projectile  be  confined  to 
the  interior  of  a hollow  cylindrical  chamber.  In  these  cases,  it  would 
be  convenient  to  mount  the  film  on  the  inner  cylinder  wall  with  the  radi- 
ation entering  the  chamber  through  a relatively  small  aperture. 


It  would  be  desirable  to  have  a single  simple  means  of  converting  a 
pair  of  independent,  simultaneous  projections  of  an  object  to  its  three- 
dimensional  coordinates,  regardless  of  the  source  and  film  configuration. 
This  report  addresses  a practical  solution  to  this  problem.  Although 
the  basic  approach  to  be  employed  is  not  original  to  this  author,2  it  is 
felt  that  the  universality  and  simplicity  of  the  analysis  developed  make 
it  a useful  tool  for  experimentors  needing  spatial  coordinates  of  pro- 
jectiles in  confined  ranges.  In  addition  to  the  presentation  of  the 
general  analysis,  a working  computer  code  and  straightforward  applica- 
tions to  two  sample  problems  are  presented. 


Superscript  numerals  refer  to  references  found  on  page  21. 
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Figure  1.  Schematic  for  Orthogonal  Projections 
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II.  THE  GENERALIZED  APPROACH 


In  the  generalized  approach  no  particular  constraints  will  be 
placed  on  either  the  location  of  the  radiation  sources  or  the  location 
of  the  film  images  as  long  as  all  locations  can  be  given  with  respect 
to  the  same  coordinate  system.  Figure  2 depicts  the  geometry  for  two 
independent  projections  of  a projectile  where  SI  and  S2  are  the  radi- 
ation sources  and  FI  and  F2  are  the  corresponding  film  images.  The 
Cartesian  components  of  these  locations  are  denoted  by  the  subscripts 
i and  j with  the  interpretation  i,j  = x,y,z.  In  this  report  vectors 
are  denoted  by  underlining,  and  the  convention  employed  in  this  analysis 
is  to  sum  over  repeated  indices  which  appear  in  a single  term.  For 
example, 

51.  F2 . = SI  F2  + SI  F2  + SI  F2  , 

11  xx  yy  zz 

and 

52.  F2.  FI.  = (S2  F2  + S2  F2  + S2,F2J  FI.. 

JJi  xx  yy  \ 8/  i 

For  convenience  we  introduce  the  unit  vector  along  the  direction 
from  j31  to  _F1  as  nl_,  and  similarly  the  unit  vector  from  to  F2  is 
denoted  by  n2^.  In  terms  of  Cartesian  components,  nl_  and  n2  are  given 
by 

nli  = CFa_-si ) . [ (Fi - sy . (Fi  -si ) j ]'1/2 

n2i  = CF2-S2J)±  [ CF2.-S2P ^ (F2_-S2J)^  ] ~1/2  . 

The  location  of  the  projectile  (i.e.,  a distinguishable  point  on  the 
projectile)  is  the  point  where  the  lines  along  nl_  and  n2_  intersect. 

One  would  find,  however,  that  the  slightest  error  in  any  of  the  measure- 
ments would  most  likely  result  in  the  lines  not  intersecting.  Con- 
sequently, the  location  of  the  projectile  will  be  obtained  by  finding 
the  place  where  the  two  lines  along  nl_  and  n2^  are  closest  to  each  other. 

The  precise  procedure  is  as  follows.  Starting  at  the  point  Sl_ 
let  us  move  along  the  direction  nl_  a distance  Dl,  and  similarly  from 
we  move  along  n2_ a distance  D2.  The  desired  location  is  given  by 
the  values  of  the  parameters  Dl  and  D2  which  minimize  the  distance 
between  the  points  on  the  lines  nl_  and  n£  corresponding  to  the  distances 
Dl  and  D2.  The  convention  for  defining  the  projectile  location  is  to 
take  the  midpoint  of  the  line  connecting  these  two  closest  points  on 
nl  and  n2.  This  convention  is  illustrated  in  Figure  3. 
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Figure  2.  The  Geometry  of  the  Generalized  Approach 


S2 


KO 


Figure  3.  The  Convention  for  Defining  the  Projectile  Location 


III.  THE  CALCULATION 


The  Cartesian  components  of  the  two  closest  points  on  nl^  and  n2 
are  given  by  : 


Dl.  = Dl  • nl.  + SI.  (1) 

111 

and 

D2.  = D2  • n2.  + S2.,  (2) 

ill 

where  the  subscripted  variables  are  Cartesian  coordinates  and  the  non- 
subscripted  variables  are  invariant  scalars.  In  order  to  find  the 
distance  of  closest  approach  we  shall  minimize  the  distance  between 
Dl  and  1X2  with  respect  to  the  independent  parameters  Dl  and  D2 . The 
distance  to  be  minimized  is  denoted  by  M: 

M <=  [(Dl*nl.+Sl.  - D2*n2.  - S2. ) (Dl*nl . +S1 . - D2*n2.  - S2.)]1/2 

In  order  to  minimize  M we  have  to  differentiate  with  respect  to  Dl  and 
D2.  In  doing  this,  one  finds  that  these  derivatives  of  M are  propor- 
tional to  1/M,  and  hence  will  be  indeterminate  if  the  lines  do  in  fact 
intersect.  To  avoid  testing  for  this  indeterminacy  it  is  convenient 
to  minimize  the  square  of  the  distance  between  Dl_  and  D2^.  For  com- 
pleteness M2  is  given  by: 

M2  * Dl2  + D22  + 2*D1  (Sl-S2;i  n^  + 2*D2  (S2-ST)i  n2i 
- 2*Dl*D2*nlin2i  + (S^-S^fSl-S^ 


The  test  to  be  used,  therefore,  is 


3M2 

3D1  302 


0. 


Since  Dl  and  D2  are  independent,  we  have  the  two  conditions 


Explicitly, 


and 


2 

0 = = 2 • D1  - 2*  D2  • nl.n2.  - 2SL.nl. 

3D1  X l l l 


jm2 

o = = 2 • D2  - 2*  D1  *nl.n2.  + 2SL.n2. 

3D2  l i xi 


where  SL  = Sl_  - S2_  . 

Solving  equations  (3)  and  (4)  for  D1  and  D2  gives 


D2  * D1  *nl.n2.  + SL.n2. 

xx  xx 


in  terms  of  Dl,  and 

D1  = 


(SL.n2.)(nl.n2.)  - SL.nl. 
v x x^ v j 3 x x 

1 - (nl.n2.)2 
X X 


(3) 

(4) 


It  is  noted  that  nl.n2^  is  simply  the  cosine  of  the  angle  between  the 
two  rays  lU  and  n^.  Therefore,  by  defining  the  sums 

COS  = nl.n2. 

x x 

SL1  SL.nl. 

x x 

and  SL2  = SL.n2.  , 

x x 

Dl  and  D2  are  given  by 

. (COS) (SL2)  - SL1 

1 ‘ 1 - (COS)2  (5) 

and 

D2  = (Dl) (COS)  + SL2  . (6) 

Since  3 (Dl)/3 (COS)-*-  ® as  COS+  1,  and  3(D1)/3(C0S)  is  a minimum  for 
COS  = 0,  the  determination  of  Dl  and  D2  will  be  least  sensitive  to 
measurement  errors  when  nl_  and  n2_  are  orthogonal . 

Substituting  the  values  of  Dl  and  D2  obtained  from  equations  (5) 
and  (6)  into  equations  (1)  and  (2)  gives  Dl_  and  D2^.  The  remaining  step 
in  the  calculation  is  to  obtain  the  coordinates  of  the  midpoint,  P,  of 
the  line  segment  connecting  Dl_  and  D2_.  The  coordinates  of  P_  are  obtained 
by  arithmetically  averaging  the  corresponding  coordinates  of  Dl^  and  D2: 

P.  = (1/2)  (Dl.  + D2±)  . 
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Using  equations  (1)  and  (2)  we  have 

P.  *(1/2) (Dl*nl.  + D2*n2.  + SI.  + S2.) 

1 X 1 1 i' 

which  is  the  desired  result. 


IV.  SAMPLE  APPLICATIONS 

As  the  analysis  was  developed  we  saw  that  the  only  input  data 

were  the  coordinates  of  the  radiation  sources  and  the  coordinates  of  the 

• * » ■ • 

film  images.  Since  the  radiation  sources  are  fixed  in  space,  their 
locations  are  constant  and  need  to  be  determined  only  once  for  a given 
experimental  setup.  The  film  image  locations,  however,  must  be  deter- 
mined from  two-dimensional  pictures  which  have  been  removed  from  the 
apparatus.  This  is  easily  accomplished  by  having  a fiducial  mark  on 
the  film  and  an  analytic  relation  between  distances  in  the  film  plane 
and  the  three-dimensional  coordinates  in  the  range  setup.  The  fiducial 
marks  will  generally  be  made  by  opaque  objects  such  as  wires  or  metal 
tabs  fixed  to  the  aparatus  at  the  exposed  surface  of  the  film. 

The  first  and  simplest  example  is  the  case  of  orthogonal  film 
plates.  For  the  particular  choice  of  Cartesian  axes  and  orthogonal 
fiducial  lines  shown  in  Figure  1,  the  translation  from  film  images  to 
spatial  coordinates  is  straightforward.  We  shall  consider  Film  1 only 
in  order  to  illustrate  the  procedure.  In  practice  one  measures  the 
distances  from  the  image  to  the  fiducial  lines.  The  x-coordinate  is 
the  perpendicular  distance  from  the  z -directed  fiducial  plus  the  distance 
from  the  z-directed  fiducial  to  the  z-axis.  Similarly,  the  z-coordinate 
is  the  perpendicular  distance  from  x-directed  fiducial  plus  the  distance 
between  the  x-directed  fiducial  and  the  x-axis.  Lastly,  for  any  point 
on  Film  1 in  Figure  1 the  y-coordinate  is  the  constant  zero. 

A simple  variation  on  the  above  example  is  the  case  where  physical 
constraints  do  not  allow  one  to  set  up  the  film  planes  at  right  angles 
to  each  other.  As  illustrated  in  Figure  4 the  coordinate  system  can  be 
chosen  such  that  one  of  the  film  planes  lies  in  a convenient  Cartesian 
plane.  Consequently,  the  x-z-plane  images  can  be  located  as  above.  To 
locate  the  image  F,  in  the  other  plane  we  shall  use  the  fiducial  lines 
labeled  1 and  2 in  Figure  4.  The  film  reading  measurements  are  the 
orthogonal  components  of  F in  the  plane  of  the  film,  F^  and  as  shown. 

The  coordinates  of  F relative  to  the  Cartesian  axes  are 
Fx  = (L1  + V C0S  6’ 

Fy  = (Lx  + Fx)  sin  0 (7) 
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Figure  4.  Nonorthogonal  Film  Planes 


and 

Fz  ■ 'W- 

where  the  convention  in  reading  the  film  has  been  chosen  such  that 
Fj  and  are  negative.  For  the  conventions  used,  the  Cartesian 

components  of  F which  are  generated  by  equations  (7)  are  used  in 
Equation  (1)  with  no  further  adaptations  required. 

The  second  example  application  to  be  presented  is  the  case  of 
cylindrical  symmetry.  As  before  the  only  requirement  is  to  relate 
the  two-dimensional  film  measurements  to  three-dimensional  coordinates 
in  the  range.  The  arbitrary  choice  for  the  Cartesian  frame  is  shown 
in  Figure  5(a). 


In  this  example  the  film  is  mounted  around  the  inner  surface  of 
a circular  hollow  cylinder,  with  inner  radius  R.  Using  the  length 
and  sign  convention  defined  in  Figure  5(b),  we  get  for  the  spatial 
coordinates  of  the  film  image  F, 


Fx  = R(!  " cos 

F = -R  sin  <J> 

y 


(8) 


where 


<J>  = _2_  radians. 
R 


Obviously  a second  simultaneous  source  and  film  image  pair  are  needed, 
and  the  same  procedure  can  be  used.  Equations  (8)  are  the  only  adap- 
tation required  in  order  to  apply  the  generalized  approach  to  the  set- 
up depicted  in  Figure  5. 


V.  THE  COMPUTER  PROGRAM 

The  computer  program  consists  of  two  subroutines.  The  first, 
SUBROUTINE  TUBES (N),  reads  and  stores  the  radiation  source  coordinates 
for  each  of  N pairs  of  discharge  tubes.  Each  pair  of  tubes  which  are 
discharged  simultaneously  is  referred  to  as  a channel.  The  sample 
programs  listed  in  Figures  6 and  7 and  the  Appendix  can  accomodate  up 
to  twelve  channels.  In  addition,  SUBROUTINE  TUBES  generates  the  com- 
ponents of  the  relative  vector  between  sources,  SL,  for  each  channel. 

The  input  for  this  subroutine  is  the  set  of  tube  coordinates  and  channel 
numbers  formatted  by  FORTRAN  Statement  101  in  Figure  6.  The  output  is 
a printed  list  of  the  source  coordinates  and  an  internal  storage  of  the 
required  vectorsaand  coordinates  in  common  block/T/. 
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Figure  5.  Film  Reading  In  Cylindrical  Symmetry 
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SUBROUTINE  TUBES (N) 

cccccccccccccccccccc 

C THIS  SUBROUTINE  READS  THE  COORDINATES  OF  THE  X-RAY  SOURCES  FOR  N 

C CHANNELS  OF  ORTHOGONAL  PAIRS  AND  STORES  THEM  IN  ISEE  FORMAT  101) 

c COHMON/T/T  112, 2, 3), S LI 12 ,3) 

8 FOR  JI  hJfR/lNNEL  NUMBER 

C FOR  J 

C 1*H0R1Z0NTAL  TUBES 

C 2*VERT IC AL  TUBES 

C FOR  K - 1*X,  2*Y,  3*Z 

(NOTF.  BOTH  THE  ORIGIN  AND  UNITS  MUST  BE  THE  SAME  AS  THOSE 
C USED  FOR  THE  FILM  READ 1NG ) «***«S L IS  DEFINED  IN  STATEMENT  1«*« 

C 

CCCCCCCCCCCCCCCCCCCC 
NRITEi  6,100) 

DO  1 K*1  ,N 

READ(5,101)NC,((T(NC,l.JI,J*1.3) ,1*1 ,2) 

NRITEI  6,101  INC,  KTINC.I  , J I , J*1 , 31  ,1  *1  ,2  ) 

DO  1 L *1  ,3 

1 SL(NC,L)*T(NC,1,L )-T (NC,2,L) 

100  FORHAT (1H1 ,12X  »17H  HORIZONTAL  TU BES , 25X , 16HV ERT I CAL  TUBES, /,8H  CHA 
CNN  EL  «7X,1HX,9X  ,1 HY ,9X ,1HZ ,19X ,1HX ,9X ,1HY ,9X , 1HZ, / ) 

101  F0RMAT(4X,12,4X,3F10.4,10X.3F10.4) 

RETURN 

END 


Figure  6.  A FORTRAN  Listing  for  Subroutine  TUBES 
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cccccc 

c 
c 


INPUT  (SEE  FORMAT  TOO! 

sfismwTOi'eHANNEi  NC  HAS  DISCHARGED 
FII.JI-F1LM  IMAGE  COORDINATES, 


IMAGE 
WHERE  1*1  FOR 


1-2 

J-l 


FOR 

FOR 


HORIZONTAL  TUDE 
TUBE 

Y,  3 FOR 


VERTICAL 
X,  2 FOR 


C 

c 

c 
c 

c 

cccccccccccccccccccc 

C0HM0N/X/SP0T(12«4)iXH!SS(12)/T/T(12,2,3),SL(12,3) 
DIMENSION  F(2,3),V(2,3IfC(2) 

READ  (5.100  )NC,5P0T(NC,4I , (IF (1,J),J*1,3>, 1-1,2) 

100  FORMAT (4X , 12 ,3X,7F10 .4) 

C CALCULATE  THE  UNIT  VECTORS  FROM  TUBES  TO  FILM  IMAGES 

C AT  STATEMENT  2 
CCII-O. 

CI2I-0. 

DO  1 K-1,2 


I 


C 

c 


33 


C 

C 

CCCCC 

CCCCC 

CCCCC 

CCCCC 

CCCCC 


rK,I)-TCNCtK,I> 

EliUSiSf  (Zi  ii  l1** ,K,,> 

C( 2) -SORT (C (2) ) 

00  2 K-1,2 

DO  2 1-1,3 
V(K,I)-VIK.l  )/C(K> 

CALCULATE  THE  DOT  PRODUCTS  OF  THE  UNIT  VECTORS  WITH  EACH  OTHER 
AND  WITH  SL 
DOT-XM ISS INC I-O. 

SL1-SL2-0. 

DO  33  1-1,3 
DOT-DOT+V ( 1,1 I *V (2,1  I 
SL1-SL1+V(1,I)*SL(NC,1) 

SL2»SL2«V(2,I l*SL(NC,I I 
CONTINUE 

D 1 -( SL 2*D0T— SL 1) /( 1 .-DOT *OOT  ) 

D2-D 1*D0T*SL2 

CALCULATE  THE  MIDPOINT  OF  THE  HISS  DISTANCE  VECTOR 

-THE  PROJECTILE  COORDINATES 
OUTPUT  GOES  TO  COHMON/X/SPOT ( 12,4) , 


WHERE  SPOTINC, 41-TIME 

SPOT  INC, 1I-PR0JECTILE-X 

SP0T(NC,2I-PR0JECT1LE—Y 
5PDTINC.3I-PR0JECTI LE— Z 

X?-f  (Ac  • I •!  I*01*VI1,I 1 

XMIS5(NC)«XM1SS(NCI-»(A1-A21*(A1— A2) 

XM 1 SS ( NC I -SORT (XHISS(NC) I 

RETURN 

END 


Figure  7.  A FORTRAN  Listing  for  Subroutine  FINDX 
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The  second  subroutine,  FINDX(NC),  reads  the  film  image  coordinates 
for  a channel  and  calculates  the  projectile  coordinates  at  the  time  the 
channel  is  discharged.  The  inputs  for  this  subroutine  are  common  block 
I’ll  and  data  cards  for  the  film  image  coordinates  which  are  formatted 
according  to  FORTRAN  Statement  100  in  Figure  7.  Care  must  be  taken  to 
make  sure  that  the  first  index  of  F(I,J)  in  FINDX  and  the  second  index 
of  T(I,J,K)  in  TUBES  correspond  to  the  same  source-image  pair.  It  is 
noted  that  the  names  "horizontal"  and  "vertical"  which  are  used  in 
Figures  6 and  7 simply  serve  to  identify  the  independent  source-image 
pairs  and  do  not  presuppose  any  particular  experimental  setup.  In  addi- 
tion, the  user  will  probably  want  to  include  a subroutine  to  convert  the 
film  image  measurements  to  three-dimensional  coordinates  as  discussed  in 
Section  IV. 

Subroutine  FINDX  performs  the  calculations  derived  in  Section  III 
and  stores  the  projectile  coordinates  in  the  array  SPOT.  For  many 
applications  it  would  be  convenient  to  add  a third  index  to  this  array 
in  order  to  identify  a round  number  or  a particular  fragment.  In  cases 
where  fragment  identification  on  the  two  films  may  be  ambiguous,  as  in 
reference  3,  the  miss  distance  M given  in  equation  (3)  can  serve  as  a 
discriminator  for  various  combinations  of  film  images.  For  convenience 
M is  saved  in  C0MM0N/X/XMISS(12) . 

In  their  present  forms  the  calling  sequence  for  the  subroutines 
TUBES  and  FINDX  is  as  follows.  TUBES  must  be  the  first  subroutine 
called  and  need  only  be  called  once  for  a given  experimental  setup; 

FINDX  is  called  for  each  pair  of  film  images  to  be  analyzed.  The  calling 
program  must  declare  the  common  block/X/SP0T(12,4) ,XMISS(12) , in  which 
the  results  are  stored.  Figure  8 is  a sample  calling  routine  for  five 
channels  with  input  and  output.  An  appendix  is  included  which  lists  the 
subroutines  TUBES  and  FINDX  in  BASIC  for  a Wang  2200  computer. 


VI.  SUMMARY 

In  this  report  a simple  and  universal  analysis  for  converting 
independent  stereo-photographic  projections  of  a projectile  to  its 
three-dimensional  coordinates  has  been  developed.  The  input  require- 
ments are  the  radiation  source  and  film  image  coordinates.  Additionally, 
sample  computer  programs  for  carrying  out  the  calculations  have  been 
presented.  These  programs  can  be  easily  adapted  to  obtain  projectile 
axes  by  locating  two  points,  and  similarly,  planes  can  be  obtained  by 
locating  three  points. 
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vi  * u*  vn 


COMMON /X /SPaT!12,6)*XHISS!12) 

READ  !5  *101  IN 
CALL  TUBES (N  I 
WRITE!  6*1001 
00  1 1*1 *N 
CALL  FINDX(NC) 

WRITE!  6,101  INC, 1SP0T !NC,J1»J*1,41*XH1SS!NC) 

1 CONTINUE 

100  FORMAT  ! 1 HI  ,8H  CHANNEL.7X  ,1HX  ,9X , 1HY, 9X ,1  HZ ,6X,6HTIHE  *6X,4HKISS) 

101  FORMAT  !4X  , 12  ,4X,6F  10  «6,E1 1 *2  ) 

END 


(a) 


0. 

500. 

0. 

300  . 

200. 

0 

100. 

500. 

200. 

300. 

200. 

100 

100. 

500. 

200. 

300  . 

200. 

200 

200. 

500. 

600. 

300  . 

200. 

300 

200. 

500. 

600. 

300. 

200. 

600 

1 

98  • 

71.4 

0. 

35.7 

•300. 

80.0 

59.9 

2 

174. 

92.3 

0. 

46.2 

•300. 

145.4 

100.3 

3 

248  • 

-66.6 

0. 

156.3 

•300. 

200.1 

149.9 

4 

353. 

118.2 

0. 

127.2 

•300. 

246.1 

207.7 

5 

498. 

•299.8 

0. 

250.0 

•300. 

285.9 

271.3 

(b) 

HOR I 70N JA L TUBES 

vertical 

TUBES 

CHANNEL  X 

Y 

z 

X 

Y 

Z 

1 

o.oooo 

500 .0000 

o.oooo 

300. 0000 

200. 0000 

0,0000 

? 

-ico.oooo 

500.0000 

200,0000 

3oo.nnoo 

200.0000 

loo.oooo 

3 

1 no .cnoo 

500.0000 

200.0000 

300.0000 

200.0000 

200,0000 

4 

-20c.onoo 

500.0000 

400.0000 

300. noon 

200.0000 

300,0000 

5 

200.onoo 

500.0000 

400.0000 

300.0000 

200.0000 

400. 0000 

CHANNEL 

X 

Y 

Z 

TIME 

MISS 

1 

49.97 

149.99 

24.97 

98 

2.97E-02 

2 

24.99 

174.99 

100.08 

174 

1.10E-ni 

3 

0.05 

200.04 

174.96 

248 

2.92E-02 

4 -24.99 

224.98 

249.98 

353 

5.38E-02 

5 

•49.81 

250.07 

324.99 

498 

5.89E-02 

(c) 

Figure  8.  (a) 

Calling  Program,  (b) 

Input  List,  (c)  Output 
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APPENDIX 


This  appendix  lists  the  subroutines  TUBES  and  FINDX  and  a calling 
program  for  these  subroutines  in  BASIC.  This  version  of  the  program  is 
operational  on  a Wang  2200  computer  with  the  input  instructions  and  out- 
put shown  on  the  CRT  screen.  The  pertinent  changes  in  notation  are  that 
SP0T(I,J)  and  XMISS(I)  are  given  by  P(I,J)  and  M(I),  respectively. 


1 REM— CALLING  AND  OUTPUT  PROGRAM— 

10  COM  T(12,6),S(12,3),P(12,4),M(12):PRINT  HEX  (03) 

20  INPUT  "NO.  OF  CHANNELS", N 
30  G0SUB  ' 100 (N) 

40  FOR  C=lTO  N : GOSUB  '31(C) :NEXT  C 
50  PRINT  HEX (03) : PRINTUSING  60:PRINT 

60  %CHANNEL  X Y Z TIME  MISS 

70  % ##  -#####.##  -####.##  -####.##  -####.##  -#.##++++ 

80  FOR  C=lTO  N: PRINTUSING  70,C,P(C,1) ,P(,2) ,P(C,3) ,P(C,4) ,M(C) :NEXT  C 
90  GOTO  330 


100  REM— SUBROUTINE  TUBES— FOR  N CHANNELS 
110  DEFFN '1 00 (N): SELECT  PRING  005:PRINT  HEX(03) 

120  FOR  K=1T0  N: PRINT  "COORDINATES  FOR  CHANNEL" ; K: INPUT  "X1,Y1,Z1,X2,Y2,Z2" 
T(K,2) ,T(K,3) ,T(K,4),T(K,5) ,T(K,6) 

130  FOR  L=1T0  3:S(K,L)=T(K,L) -T(K,L+3) :NEXT  L 
140  NEXT  K: RETURN 


200  REM SOUBROUTINE  FINDX FOR  CHANNEL  C 

210  DEFFN' 31 (C) : SELECT  PRINT  005:PRINT  HEX(03):DIM  F(2,3) ,V(2,3) ,C(2) 

220  PRINT  "FOR  CHANNEL" ;C:PRINT 

230  INPUT  "X,Y, Z FOR  FILM  IMAGE  1",  F(1 , 1) ,F(1 ,2) ,F(1, 3) : INPUT"X, Y,Z 
FOR  FILM  IMAGE  fi",F(2,l),F(2,2) ,F(2,3) :INPUT  "TIME",P(C,4) 

240  C(L)  ,C(2(»0 

250  FOR  K=1T0  2:FOR  I=1T0  3:V(K,I)=F(K,I)-T(C,I+3*(K-1)) :C(K)=C(K)+V(K, I) 
+2  :NEXT  I -.NEXT  K 

260  C(1)=SQR(C(1))"C(2)+SQR(C(2)) 

270  FOR  K=1T0  2:FOR  I=1T0  3:V(K, I)/C(K) :NEXT  I:NEXT  K 
280  D,S1,S2,M(C)=0 

290  FOR  I»1T0  3:D=D+V(1,I)*V(2,I) :S1=S1+V(L,I)*S(C,I) :S2=S2+V(2 , I) *S (C, I) : 
NEXT  I 
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300  D1=(S2*D-S1)/(1-DI2) :D2-D1*D+S2 

310  FOR  I-1TO  3 : A1=T(C , I) +D1*V (1,1) : A2=T(C , 1+3) +D2*V (2,1) : P (C , I)  = (A1+A2) / 2 
= M(C)+(A1-A2)+2:NEXT  I 
320  M (C)  =SQR  (M (C)  ) : RETURN 
330  REM 
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